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We study excitations of the predicted antiferromagnetically ordered vortex cores in the supercon- 
ducting phase of the newly proposed SO(5) model of strongly correlated electrons. Using experi- 
mental data from the literature we show that the susceptibilities in the spin sector and the charge 
sector are nearly equal, and likewise for the stiffnesses. In the case of strict equality SO(5) symme- 
try is possible, and we find that if present the vortices give rise to an enhanced neutron scattering 
cross section near the so called tt resonance at 41 meV. In the case of broken SO(5) symmetry two 
effects are predicted. Bound excitations can exist in the vortex cores with "high" excitation energies 
slightly below 41 meV, and the massless Goldstone modes corresponding to the antiferromagnetic 
ordering of the core can acquire a mass and show up as core excitation with "low" excitation energies 
around 2 meV. 



I. INTRODUCTION 

Inspired by the discovery of a sharp antiferromagnetic 
resonance, later denoted the tt resonance, in neutron 
scattering experiments on the superconductinffi-chase of 
YBa2Cu3C>7 at (tt, it) in the reciprocal spaceEra a new 
idea was introduced a year ago in the search for a theory 
of the high-T c superconductors. Within -the framework 
of the t-J limit of the Hubbard model&Q it was shown 
theoretically already three years agotl that the tt reso- 
nance could be explained in terms of a new collective 
mode in the particle particle channel of the model, and 
that this mode is intimately connected with the symme- 
try of the superconducting gap. Pursuing the symmetry 
aspects of the problem Zhang proposed a theory combin- 
ing antiferromagnetism and superconductivity by sym- 
metry arguments.tj The operators responsible for the tt 
resonance was identified with the six generators of ro- 
tation between the antiferromagnetic state and the su- 
perconducting state. Furthermore, it was proposed that 
the phase diagram of the cuprates can be understood 
as a competition at low temperatures between d-wave 
superconductivity and antiferromagnetism of a system 
which at higher temperatures posses SO (5) symmetry. 
The group 50(5) is sufficiently large to accommodate 
both the gauge group U(l) (= 5*0(2)) which is broken in 
the superconducting state, and the spin rotation group 
5*0(3) which is broken in the antiferromagnetic state. In 
the simplest version the transition between the two or- 
dered states is controlled by one parameter - the chemical 
potential for holes in the otherwise half-filled quadratic 
lattice of spin 1/2 fermions. Tie-idea has generated a lot 
of discussion among theoristsoil3, and no consensus on 
the matter has emerged. 

We shall not enter this discussion here. Rather we 
will take the approach of assuming the model to be a 
fair description of the high-T c materials and derive ex- 
perimental consequences, which can be tested in the lab- 
oratory. Arovas et al.c3 have pointed out that in the 



vortex cores of fluxoids in the superconducting state the 
order parameter will escape into the antiferromagnetic 
subspace, meaning that in these cores we have local an- 
tiferromagnetic moments instead of a simple featureless 
normal metal core. This unique prediction of the 50(5) 
model should in principle be quite simple to verify exper- 
imentally. However, preliminary measurements^] where 
one looked for Bragg-scattering from these moments be- 
longing to the vortex cores did not produce any signal. 
This is perhaps not so surprising, since each of the vor- 
tex cores will form a one-dimensional magnet (along the 
c axis) and at finite temperatures such a system does 
not form long range order, and no Bragg peak is to be 
expected. 

In this paper we are going to pursue the idea that in 
each copper-oxide plane small islands of antiferromag- 
netically ordered spins exist associated with the vortices 
generated by an external magnetic field. The direction 
of the spins in these islands will not be very strongly cor- 
related from layer to layer and from island to island. In 
one island, however, there should exist excitations of the 
spins, a kind of bound spin wave modes or size quantized 
magnons. Using samples in the under-doped regime of 
the superconducting phase, where the proximity of the 
antifcrromae*S|etic phase stabilizes the antiferromagnetic 
vortex corestil, one should be able to pick up these core 
excitations in inelastic neutron scattering measurements. 
The modes can be classified according to the approximate 
symmetry: Two zero-energy or "low" energy Goldstone 
modes whose existence is guarantied by the exact spin 
rotation symmetry, and two resonances or weakly bound 
"high" energy modes related to the tt resonance arising 
form the approximate 50(5) symmetry allowing for ro- 
tations between the d wave superconducting phase and 
the antiferromagnetic phase. Based on experimental data 
taken from the literature we discuss in Sec. [n]the 50(5) 
model and its parameters. The values of the suscepti- 
bilities and the stiffnesses in the charge sector and the 
spin sector are found to be nearly equal, a remarkable 



1 



fact supporting the 50(5) model. In Sec. Ill we set up 
the calculation for excitations of the vortex core. In the 
isotropic case we show analytically that the Goldstone 
modes of the vortex core indeed have zero energy, and 
that the vortex indeed generates a resonance reminiscent 
of the 7r resonance at the bottom of the continuum. In 
the anisotropic case the Goldstone modes remains mass- 
less, but for certain anisotropies the it resonance can be 
transformed into a bound state localized at the vortex. 
Treating the external fields more accurately introduces 
symmetry-breaking terms in the Hamiltonian and results 
in finite masses to the Goldstone modes. This will be dis- 
cussed in Sec. |Tv| . Concluding remarks are contained in 
Sec. 0. 



II. THE SO(S) MODEL AND ITS PARAMETERS 

In the £0(5) model the relevant order parameter is a 
real vector n in a five dimensional superspin space with 
a length which is fixed at low temperatures, 



n = {ni,n 2 ,n 3 ,n 4 ,n 5 }, 



= 1. 



(1) 



This order parameter is related to the complex super- 
conductor order parameter, ip, and the antiferromagnetic 
order order parameter, m, in each copper-oxide plane as 
follows: 



tp = f e ie = m + 



(n 2 ,n 3 ,n 4 ) 



(2) 



In Ref. ^ Zhang argued how in terms of the five dimen- 
sional superspin space one can construct an effective La- 
grangian £(n) describing the low energy physics of the 
t-J model. In the absence of external electromagnetic 
fields it takes the form 
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where the indices run from 1 through 5. 

The generalized susceptibilities, x, fall in three groups: 
Xc = Xi5 connecting the charge sector {ni,n^\ with it- 
self, Xs = X23 = X24 = X34 connecting the spin sector 
{n 2 , n 3 , n 4 } with itself, and Xtt = Xi(2,3,4) = X(2, 3,4)5 con- 
necting the spin sector with the charge sector. Similarly 
with the stiffnesses: p c = pis, p s = p 23 — P24 = P34, 
and p n = /Ci(2,3.4) = P(2,3.4)5- Below, based on experi- 
mental data, we find that p s ss p c and Xs ~ Xc- It is a 
remarkable fact that the dynamics in two such distinct 
sectors as the charge and the spin sector are governed by 
coupling strengths of the same size, and it can be taken 
as one strong indication of the near 50(5) symmetry of 
the cuprates. At present the values of x* and p-n are not 



known experimentally, and it is part of our work to es- 
tablish a method to measure them. If the corresponding 
coupling strengths are the same in all sectors we denote 
it the isotropic case, otherwise the anisotropic case. 

In the following we estimate on a 25 % accuracy level 
the typical zero temperature values of the parameters 
of the SO (5) model obtained for various cuprates of 
the form YBa 2 Cu 3 6+a: (YBCO) and La 2 -xSr x Cu04 
(LSCO) with a range of doping levels x. All numerical 
values are listed in Table Q. FiisJj we note that for both 
materials the Cu-O-Cu distanced is a = 3.8 A. However, 
YBCO contains two CuO planes over a distance of 11.4 A 
(the sum of the alternating inter-layer distances 3.2 A 
and 8.2 A), while LSCO contains one CuO plane every 
6.6 A. Hence, when needed for normalization purposes 
we employ c = 6.1 A as the typical inter-layer distance. 

From the upper critical field H^Xhe typical correlation 
length is found to be £ ~ 16 AJlJ while many different 
methods like muon spin rotation, magnetic torque, mag- 
netization, and kinetic inductance (see Ref. [l5| and refer- 
ences therein) all result in a London length, Xl — 1350 A. 
Thus the Ginzburg-Landau parameter is k = A^/£ ~ 84. 

The connection between the generalized coefficients 
and the directly measurable parameters are given below. 
In the completely isotropic case where all generalized co- 
efficients are equal we have:Eil 



A L = Y e V c /ppo- 



(4) 
(5) 

(6) 



To estimate p s and Xs of the spin sector, experimental 
measurements are combined with theoretical calculations 
of spin waves within the two dimensional spin li2_quan- 
tum Heisenberg model of antiferromagnetism.Ejlij The 
bare coupling constant J is related to p s , Xs and the spin 
wave velocity v s asc£l 

Ps = Z p J/A (7) 
Xs = Z x p 2 a /8J (8) 

V s = Z c y/p/x = ZcVZJ/Pa, (9) 

Z p = 0.72, Z x = 0.51, Z c = 1.18, (10) 

where for brevity a momentum p a = ti/a has been in- 
troduced, and where the Z's are renormalization con- 
stants, which for classical spin waves all equals 1, but 
differs from 1 when quantum fluctuations and spin wave 
interactions are taken into— account. Neutron scatter- 
ing experiments on LSCOlJ have led to a determina- 
tion of v s and from that to J = 132 meV in agree- 
ment with other experiments. From J one calculates 
p s = 24 meV. Independently, p s have been determined 
by neutron scattering measurementE3 of the antiferro- 
magnetic correlation length £,afm(T) cx exp(2Tr p s /ksT) 
also leading to p s — 24 meV. From Eq. (0) one finds 
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-53 peV/A 2 


ps ^ 


0.024 eV 


Xs/pl ^ 


0.49 eV -1 


5 P > 


-1 


Pc ^ 


0.018 cV 


0.42 eV -1 


<5 X > 


-0.014 



TABLE I. The average values of typical parameters of 
the 50(5) model based on experimental data for YBCO and 
LSCO as described in the text. The Cu-O-Cu distance is 
denoted a, while the average distance between the CuO planes 
is denoted c. For brevity a momentum p a = h/a has been 
introduced. The number of spins in a given vortex core is 
estimated by Tv£ 2 / a 2 . T he anisotropy parameters S p and S x 
are defined in Sec. [II B. 



Xs/pi = 0-49 eV 1 . Ap^pore recent neutron scatter- 
ing experiment on LSCOEil yielded the consistent result 
Z x = 0.4 ±0.1 and J = 125 meV. 

In the charge sector we have p c = (c/ po)(1i/2e\L) 2 , 
while Xc — Pc/v 2 , where v c is the sound velocity of the 
electron liquid. The charge stiffness is readily obtained 
from the London length, p c = (h/2e) 2 (c/ '/zo)(1/Al) 2 w 
18 meV. The susceptibility Xc is obtained as follows. 
In ordinary superconductors the Goldstone mode corre- 
sponding to the sound velocity is rendered massive by 
the Anderaan-Higgs mechanism and turned into a plas- 
mon modeE3, but as a theoretical concept it can be cal- 
culated, and it is found to bC|-afj-Lhe order of the Fermi 
velocity, vp. Detailed studiesE3a have shown that the 
dispersion relation for quasi particles moving around on 
the lattice of the t-J model at low doping near half fill- 
ing is governed not by t but is renormalized t* sa J. 
Using the dispersion relation (3.15) of Ref. [ll|, £k = 
const + Jcos(k x a) cos(k y a) + | J[cos(2k x a) + cos(2fc y a)], 
yields v c « Vp w 1.58J/p a . From this we find Xc/Pa ^ 
0.42 eV^ 1 . 

To estimate the value of the phcnomcnological symme- 
try breaking constant g in Eq. (|^) we consider complete 
antiferromagnetic ordering, i.e. |m| 2 = 1. In this case 
the ordering energy density in the SO(b) model is simply 
— |p|m| 2 = — \g. On the other hand this energy density 
can also be expressed within the t-J model as — J/(2a 2 ), 
and therefore g = J /a 2 ~ 9 meV A -2 . The effective 
coupling constant g of Eq. (Q) is much smaller. Antic- 
ipating the discussion in Sec. Ill A of the it resonance 
frequency uj w rj 41 meV we find the following estimate: 
g = - X ^l ~-53/ieV/A 2 . 



III. VORTICES AND CORE EXCITATIONS IN 
THE SO (5) MODEL 

In Ref. |ll] the vortex solutions to the isotropic SO(5) 
model have been studied in great detail. To study the 
anisotropic case we are going to use a different method. 
To est ablish our method and notation, we will start out 
in Sec III A by deriving some of the known results for the 



isotropic case, before in Sec [II B we continue with the 
anisotropic case. 



A. Vortices and resonances in the isotropic case 

In the symmetric version of the SO(5) model the gen- 
eralized susceptibilities and stiffnesses are isotropic in su- 
perspin space, and the only symmetry breaking terms are 
quadratic terms governed by the chemical potential p, and 
the phenomenological constant g chosen such that su- 
perconductivity is favored. The external electromagnetic 
fields will now by included through the vector potential 
A. However, in this section we only keep the interac- 
tion with the tp part of the order parameter this being 
the dominating part of the external fields. In Sec. IV A 
we will include the Zeeman interaction between A and 
the m part of the order parameter and demonstrate ex- 
plicitly that this only leads to minor changes. In this 
approximation the Lagrangian £j SO then has the form: 



\x\dtM 2 
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where c is the lattice constant perpendicular to the 
copper-oxide planes and c the speed of light. For later 
convenience we have incorporated the constraint Eq. (|l|) 
through the Lagrange multiplier A. Using dimensionless 
polar coordinates (s,4>), with s = r/£, centred at the 
origin of the vortex core, we seek solutions of the form 



^M)=/(*)e-*, 



m(s)ep, A 



2e~l s 



(12) 



where is an arbitrary unit vector in (712, n^, ri4)-space 
(equivalent to real space) taken to be (0, 1,0) in the fol- 
lowing, while is the azimuthal unit vector. The Euler- 
Lagrange equation for A yields the constraint 



/(s) 2 + m( S ) 2 = l, 



(13) 



such that once f(s) is determined so is m(s). The Euler- 
Lagrange equation for m is used to express A(s) in terms 
of m(s): 



1„ 2 (A- 5 )£ 2 
m p 



(14) 



Eqs. (H) and (0) are then used to eliminate m(s) and 
A(s) in the Euler-Lagrange equation for f(s), and as in 
Ref. |ll| we end up with: 
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The Euler-Lagrange equation for a(s) becomes 



<9 c a <9 s a = 



(a-1) 



/ 2 



= 0. 
(15) 

(16) 



Equations (|T^) pad (16) are solved by the numerical 
shooting method!!^ and yields f(s) and a(s). 

Introducing the "effective potential" Vq (s) as the right 
hand side of equation Eq. (HJ) , 
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(A- 5 )e 2 



— V 2 TO 
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fV 2 f 
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we can use the solution of f(s) to determine Vq(s), which 
in turn results in the Euler-Lagrange equation for m of 
the form 

[-V 2 + F (s)] 777 = 0. (18) 
Naturally, by construction, Eq. (|l8|) , is automatically ful- 



filled using m(s) = yl — /(s) 2 , but for the forthcoming 
studies of core excitations it is useful to think of the static 
core as corresponding to the zero energy solution of the 
Schrodinger-like equation Eq. ( |l8| ) where Vq(s) clearly 
plays the role of an effective potential. For the same rea- 
sons we rewrite Eq. ( p~5| ) for /(s) by the use of Eq. Jl7j ) 
forVb(s): 



-V 2 



Vo(s) 



a(2 - a) 



/ = /• 



(19) 



Thus having established the notation and found the 
static vortex solutions we now turn to the problem of 
finding excitations in the vortex core. We denote the 
static vortex solution by no, and due to the SO(3) sym- 
metry in the spin sector we are free to choose the di- 
rccti on of no at will. In anticipation of the treatment in 
Sec. EVA where an external magnetic field in the z direc- 



tion forces no to lie in the xy plane we choose no oc e,,, i.e. 
only the second of the three no components is nonzero: 

n = {/(s) cos(0), 0, m(«), 0, -/(«) sin(</>)}. (20) 

We seek excitations Sn which are of lowest order in the 
deviations Sip and Sm. These turns out to be perpendic- 
ular to no, i.e. Sn cx e x or Sn oc e z or linear combinations 
thereof, and hence of the form 



Sn x 
Sn z 

Sm 



= {0,1, 0,0,0} Sm, 
= {0,0, 0,1,0} Sm, 



Sm n i(s)e 



14 f 



(21) 



Throughout this work we are dealing with cylindrical 
symmetric vortices, so the excitations are characterized 



by the angular momentum I and the radial index n. 
Of course, Sn is not a complex vector, so the notation 
exp(il(p) is merely a short hand notation for either cos(l<p) 
or sui(l(p). The Lagrangian is now written to second or- 
der in Sn as C = £)^a( n o) + C^(Sn). The explicit form 
of the second order term is 



CW(Sn) = ± X \d t Snf 
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p\VSn[ 
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\g-\{a)] \Sn{ 



(22) 



Assuming solutions Sn of the form given in Eq. ( |2l| ) the 
Euler-Lagrange equation for 5m n i(s) then takes the form 
of the following eigenvalue equation: 



c 2 



Vo(s) 



Sm n i = eSm n i, 



(23) 
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X 2 
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Using the approximate SO(5) symmetry we can im- 
mediately find two analytical solutions <5tooo an d Smoi 
to Eq. (|23[) . Due to the exact spin rotation symme- 
try it does not cost any energy to rotate the order pa- 
rameter no of Eq. ( |20| ) in spin space. Rotating no a 
small angle S9 in the (713,714) plane produces the devi- 
ation <5noo = {0, 0, 0, S9 m(s), 0}. If the trial solution 
(5 777.00(5) tx m(s) is used (note that I = 0), we see from 
Eq. ( |is| ) that as expected Eq. fl23| ) is satisfied with e = 0. 
Thus in the the effective potential description the ground 
state vortex configuration corresponds to a zero energy 
and zero angular momentum mode. 

The second solution is found by rotating between the 
charge sector and the spin sector. In a perfect 50(5) 
symmetric model such a rotation does not cost any en- 
ergy. However, one central idea in the SO(5) model of 
high T c superconductors is that the SO(5) symmetry is 
only approximate. It costs a finite energy to rotate be- 
tween the spin and charge sectors. Experimentally this is 
reflected by the 41 meV ir resonance, and theoretically by 
the symmetry breaking term \gm? of Eq. (^). We thus 
expect that by rotating no a small angle 50 in the (77,1, 774) 
plane a deviation Sn^i — {0, 0, 0, SO f(s) cos(</>), 0} is pro- 
duced which is an eigen-excitation with e > 0. If the trial 
solution Smoi(s) oc /(s)cos(^) is inserted into Eq. 
(note that I — 1), we obtain 



1 



+ V (s) 



f = ef. 



(24) 



It is seen from Eq. (|l^) that this equation is indeed sat- 
isfied by f(s) if e = 1 and if the additional potential 
a(2 — a)/s 2 can be neglected. It turns out that the large 
value of the Ginzburg-Landau parameter, k ~ 84, indeed 
does make the additional potential negligible. Numeri- 
cal calculations show a(2 — a)/ s 2 < 0.001 (for any value 
of s) which is much smaller than Vq(s) ~ 1 and e = 1. 
The approximate nature of the solution is not surprising, 
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since the external magnetic field does in fact break the 
50(5) symmetry by coupling only to the (nijUs) com- 
ponents of the order parameter. However, the larger a 
value of k the weaker this symmetry breaking appears, 
and in the limit of infinite n the approximate solution 
becomes exact. Since the ground state energy is set to 
be zero as the zero energy mode, the excitation energy 
of the dmoi resonance is given by e = 1, or going back 
to frequency: = y/—g/Xi which is in accordance with 
that of the ir resonance given in Ref. ||. 

In conclusion we note that both <5moo and <5tooi con- 
tains no radial nodes, hence the notation n = 0. Any 
other excitation or resonance would contain more nodes 
and thus have higher energies. Since 8m 01 corresponds 
to a resonance at the bottom edge of the continuum we 
can infer that for the isotropic case no bound collective 
excitations exist in the antiferromagnetic vortex core. 



B. Vortices and excitations in the anisotropic case 



V 2 to (A - g)e 



m p n 

d^a — -d s a H 

s 



1 + 6. 



[(Vm) 2 +mV 2 m] , 



-/ 2 



1 - 



(28) 
(29) 



It is seen how the anisotropy leads to more nonlinear 
terms in the differential equations. Using the same nu- 
merical methods as in Sec. [II A we study the static vor- 
tex cores for various values of S p and 5 X . Some results 
are shown in Fig. |l|. 





We now turn to the anisotropic case. As discussed 
in Sec. || p s w p c and Xs ~ Xc- In the following we 
therefore study the consequences of anisotropies arising 
from P7T ^ p s = p c and Xtt ^ Xs = Xc 



P = Ps = Pc 

Pit = p+Ap, 

Sp = Ap/p, 



X = Xs = Xc 
X* = X + A X , (25) 



Ax/x- 



These anisotropies are not known experimentally, and it 
is part of our work to establish a method to measure it. 
We begin by finding constraints on them. The stiffness 
has to be a positive number, i.e. p„ = p(l + S p ) > 
or S p > — 1. The n susceptibility Xn is related to the 7r 
resonance and to the coupling strengths g and g. Defining 
7 = (g— go)/ g, where the subscript refer to the isotropic 
case, we can estimate 7 using the values listed in Table | 
and find 7 ~ 1.014. The relation between g and go can 
be written as 
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2P,2 = l-(l + * x )7 



S'-X7r(-r-) =5o 



1-7 



(26) 



To ensure the superconducting phase it is mandatory to 
have g < and hence from the enumerator in Eq. (|2C 



that 1 - (1 + 5 x )j < or -0.014 < S x . This constraint 



is listed in Table | 

The Lagrangian £ ani = £ iso + AC in the anisotropic 
case differs from the Lagrangian £j so in the isotropic case 
by AC containing terms proportional to the anisotropies 
8 P and S x . The Euler-Lagrange equations for /, to, and 
a corresponding to C anl become 



V 2 / 



f(dsf) 2 
I-/ 2 

1 +5v 



l-(- 



1, 



a-!) 



2S n 



/(I - / 2 ) (27) 

f(i-n 



5 P = 0.5 
S p = 0.0 
8„=-0.5 




FIG. 1. A display of the vortex cores for various values 
of the anisotropy parameters. Graphs of the modulus f(s) of 
the superconducting order parameter as well as of the effective 
potential V(s). The parameters are k = 84 and 7 = 1.014. 
When 5 P is varied, 5 X — and vice versa. 

In the isotropic case the static vortex core led to an 
effective potential description with Vo(s) defined as the 
right hand side of Eq. (n4|) . Similarly, for the anisotropic 
case we now define an effective potential V(s) as the right 
hand side of Eq. (|2§|) which then satisfies the equation 



-V 2 + V{s)} to = 0. 



(30) 



V(s) can be expressed in terms of / by using to 2 = 1 — / 2 , 
and thus it can be found by solving Eqs. ( pq ) and (|2S|). 
Due to the anisotropy the effective potential and hence 
the excitation spectrum changes. The effective potential 
for some anisotropy parameters are shown in Fi g, pi . We 
study the transverse excitations given by Eq. ( |2l|) and 
find the following eigenvalue equation to be fulfilled by 
Sm n i(s): 



r 2 



V(s)] Sm r , 



Pt, 



uj Sm n i = e 5m n i. (31) 



As in the isotropic case a zero energy solution corre- 
sponding to the static vortex solution is trivially given. 
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What is new in the anisotropic case is that now bound 
excitations do exist with e < V^oo). Some of them are 
shown in Fig. g. 




FIG. 2. The radial part 5m n i of bound antiferromagnetic 
excitations in the vortex core in an anisotropic case. To show 
the effects more clearly we have chosen the rather extreme 
parameter value 5 P — —0.5, while having k = 84, 8 X = 0, 
and 7 = 1.014. To the left is shown the zero energy state 
Smoo corresponding to the static vortex core. To the right 
are shown the first two excitations Smoi and Smio above the 
zero energy mode. The full line in both panels is the effective 
potential V(s). The ordinate axis accounts for V(s) in the en- 
ergy units of Eq. (^). The excitations 5m n i are given in arbi- 
trary units. The eigenenergies are represented by the dashed 
horizontal lines with £oo = 0.00, eoi = 0.81, and eio = 0.97. 



In Fig. |^ we plot as a function of the anisotropy pa- 
rameter 5 P the energy of the lowest excited state Smoi 
above the ever present zero energy mode <5tooo- It is seen 
how Smoi evolves from the scatt ering resonance, the it 
resonance, discuscd in Sec. Ill A with e = V^(oo) at the 
isotropic point S p = to a strongly bound state with 
e = at S p = — 1, the largest negative value allowed. At 
this point the excited state thus coincide energetically 
with the zero energy state indicating the phase transi- 
tion from superconductivity to antiferromagnetism - the 
gap of the 7r excitations has collapsed. 

Our calculations thus lead to the following predic- 
tion. For a given anisotropy S p < neutron scattering 
will in zero magnetic field show the n resonance with 
£ = Kmi(oo) or frequency u> — oj n . As the magnetic 
field is turned on more and more vortices are created. 
Each of them supports a <5moi excitation with an energy 
£oi < V^(oo) or frequency u>oi < to v . Consequently an ex- 
citation resonance should show up at the low-energy side 
of the 7r resonance, and the amplitude of this resonance 
scales with the magnetic field. 

Similar results does not hold for S x ^ 0. The reason 
for this can be traced back to the behavior of V shown 
in Fig. |l|. When S x is made positive the potential widens 
as in the case of negative S p , however, at the same time 
the asymptotic value drops. As a result no excitations are 
bound in the potential. In the case of negative S x the po- 
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FIG. 3. The energy e of the lowest bound excitation 
as a function of the anisotropy S p for k = 84, S x — 0, and 
7 = 1.014. The dotted line is a fit, e{S p ) ~ 1.0 -0.75p -0.3(5*. 
Note that e(0) = 1 and e(-l) = 0. 



tential narrows down, but although the asymptotic value 
rises, it does not rise enough to bind any excitations. 

It should be added that this change in the asymptotic 
value of V due to Sx is in accordance with the result 
for the 7r resonance in a bulk superconductor without 
the presense of vortices. In the dimensionless units of 
the problem, see Eq. (23), the frequency of the n 
resonance is given by 



Xir,,a i , 9 A X (2^) 2 
-9 9 fi 



(i- 7 ) + Vr 
(!-7) - Vy' 



(32) 



Inserting the values 5 X — —0.005, 0, and 0.005 we find 
2.14, 1.00, and 0.47, respectively, as V(oo) in Fig. 0. 



IV. GOLDSTONE MODES 



In Sec. |I| we studied the zero energy modes and the 
excitations of the antiferromagnetic order parameter m 
in the superconducting vortex core. The zero energy 
modes 5n x and <5n 2 in Eq. ( pi] ) were found to be degener- 
ate gapless Goldstone modes. This degeneracy is a result 
of two approximations, one being the neglect of the Zee- 
man interaction between m and the external magnetic 
field B, the other being the neglect of inter-layer inter- 
action between the spins. In this section we show how 
the degeneracies are lifted and massive Goldstone modes 
appear when the approximations are abandoned. 



A. Coupling of core excitations to the external 
magnetic field 

The primary effect of the external magnetic field B = 
Be z on the system is the creation of superconducting 
vortices through the interaction with the superconduct- 
ing order parameter ip. To a good approximation the 



G 



antiferromagnetic core described by m can be treated 
disregarding the Zeeman coupling between B and m. We 
now take this coupling into account. Using standard field 
theoretic methodso the spin operator Sj acting on site 
j can be expressed by a classical field hj, 



(33) 



where s = 1/2 and Q = 7r(a 1 , a 1 , c x ) is the antiferro- 
magnetic ordering vector. We then have exp[iQ • Rj] = 



(-iy 



1 on sublattice A and —1 on sublattice B. 



Since the system is close to be completely antiferromag- 
netically ordered hj is written as 



m, + (-l)'lj 



(34) 



where lj denotes a small ferromagnetic component on 
top of the antiferromagnetic background nij. Both L, 
and rrij are slowly varying fields in space. The rapid 
antiferromagnetic variation from site to site is explicitly 
taken into account by prefactors (— l) 3 . The smallest 
deviation possible is obtained by having lj perpendicular 
to nij, and since hj is normalized to unity, we obtain to 
lowest order in L that 



N 2 = i, 



Ij-mj =0, 



m., 



= 1. 



(35) 



The Hamiltonian H& corresponding to the coupling be- 
tween the spins and B is 



E 

3 



i 5 >sB-^(-l)^hj. (36) 



Only the ferromagnetic component yields a non-zero con- 
tribution to the sum, and after taking the continuum 
limit we end with a Lagrangian density Cb given by 



I^B.l(r). 



(37) 



The appearance of a non-zero ferromagnetic compo- 
nent l(r) leads to a loss of antiferromagnetic ordering en- 
ergy. This has to be included in the description through 
the Hamiltonian Hafm describing the spin-spin interac- 
tion: 



H 



AFM 



(38) 



Cm') 



{3,3') 



Performing the sum and taking the continuum limit re- 
sults in the Lagrangian density £afm, 



£afm 

The constraint L 



(39) 



nij = is incorporated through 
a Lagrange multiplier A' and we end with the following 
Lagrangian density CJ of the magnetic effects 



i j 

2^2 



l(r) 



•l(r) 



1 J 

2 a*' 



(40) 



An effective Lagrangian £ m B(m) for the interac- 
tion between m and B is found from the partition 
function Z = J X>m fVlJ VX' exp[iS(m, 1, A') /K] = 
j r Dmex-p[iS c ff(m)/Ti] from which the final result can be 
extracted after integrating out 1 and A': 



£ mB (m) = --L(B-m) 2 



-LB 2 



2^' 



(41) 



where L = (g*/is) 2 /(4Ja 2 ). Adding £ m n to the La- 
grangian of the anisotropic case £ an i leads to the follow- 
ing eigenvalue equation for an excitation Slip = Smep, 
with f3 = x or z, 



[-V 2 + V{s)]8m 



Ptt 



8m. (42) 



As was the case for Eq. ( p3[ ) two eigensolutions propor- 
tional to 77i can now be found . One mode has (3 = x 
(i.e. e z -ep — 0) and remain a zero energy mode: 



TlUJr, 



0. 



(43) 



The other mode has j3 — z (i.e. e z -e.fj — 1), resulting 
in a nonzero energy, which can be estimated by setting 
B = $o /7rA 2 and using the parameters of Table || 



TlLJ-, 




V2h 2 

m*A? 



5.9 /jeV. 



(44) 



We can thus conclude that the presense of the external 
magnetic field in fact does break the degeneracy of the 
two gapless Goldstone modes, leaving only one mode gap- 
less while rendering the other massive. The rather small 
value of the gap, 5.9 /xeV « 68 mK, would be very dif- 
ficult to observe in a neutron scattering experiment. It 
is two orders of magnitude smaller than the ir resonance, 
and hence it is seen to be a good approximation to disre- 
gard the interaction between the external magnetic field 
and the antiferromagnetic order parameter. 



B. Inter-layer coupling 

To this point we have only treated one single CuO 
plane. Naturally, to stabilize the order parameter, inter- 
action between the layers has tacitly been assumed. In 
this section we are explicitly going to include that part 
of the inter-layer coupling which arises from the antifer- 
romagnetic coupling between the CuO planes. We model 
this coupling by a Hamiltonian H' , where spins at site j 
in the CuO plane £ are interacting with the closest spins 
in the neighboring planes C ± 1 as follows: 
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(45) 



The inter-layer coupling J' is much smaller thaa the in- 
tralayer coupling J, namely J' ~ 4 x 1CP 5 J.E3 In this 
specific model there is no coupling between sites having 
different in-plane index j, and we are led to consider ID 
spin chains perpendicular to the planes. As before, see 
Eqs. ([33]) and (|34|), the spins are represented by the clas- 
sical fields h, m, and 1: 



S U — » < 1)4., 



m 



+ (-!)%■ c . (46) 



Henceforth, we drop the site index j and focus on just one 
of the spin chains. The Hamiltonian H' is now expressed 
in terms of the classical fields, and in the continuum limit 
we obtain the form 



d 2 v 



C 



{2-[h f (r)-h f _ 1 (r)] 2 }. (47) 



Defining Am^ = 1 wc obtain the following form 

of the Lagrange density, C\\ c for the inter-layer coupling: 



J' 



2E 

c 



[2 



Am?] 



(48) 



As in Sec. IV A, excitations 5m in the order parameter 
m = me y are sought in the perpendicular directions e x 
and e z . The motion in these two directions is indepen- 
dent of each another, and we write the excitations <5m^, 
(3 — x or z, as 



5m p (Q = 6f,({) e , 



(49) 



where the amplitude 9p(() is m times the (small) angle 
by which the order parameter m in plane £ is tilted away 
from its equilibrium position. The eigenmodes are found 
by the Fourier transformation 



ikC,c 



2tt 



1,2,3,. 



(50) 



where N c is the number of CuO planes in the sample. 
Because of the Fourier transform it is more natural to 
work with Lagrange functions, L — J Cd 2 r, rather than 
Lagrange densities. However, in the end by dividing 
L with N^ 2 , the number of planes times the effective 
area of a vortex, the Lagrange function is rendered into 
a Lagrange density. The resulting Lagrangian density 



CflJ(6m.p) thus becomes 

4c ) (^) = E{fi 9 ^ 

u a v 



k,0 



Aa 2 



[l-cos(frc)]|0* l/3 | 



(51) 



where, after using x 
uJk.p ars seen to be 



H 2 /8Ja 2 , the eigenfrequencies 



Uuik,p = 2y/ J'J[1 — COs(fcc)]. 



(52) 



The inter-layer coupling thus splits the 2A c -fold degener- 
ate Goldstone modes in a stack of N c vortices. The two 
k = modes, u>o,p = 0, remain zero energy modes. How- 
ever, since cJk,[x>x) — UJ 2-K-k.{x,z)i the rest of the modes, 
k > 0, split up in a quasi continuous band consisting of 
A c /2 — 1 four-fold degenerate massive Goldstone modes. 
The most massive modes are found for k = n/c with an 
energy 



TlLO 



tt/c/3 



V8J'J ~ 2.2 meV w 26 K. 



(53) 



We note that for two-layer compounds such as YBCO 
the effect of the inter-layer coupling can be estimated by 
setting N c = 2 in the above calculation. In this case only 
k = and k = n/c occurs, and the degeneracy splitting 
is given by the stated 2.2 meV or 26 K. We conclude 
that also the inter-layer coupling produces only minor 
effects in the excitation spectrum as compared to the 7r 
excitation, however, the estimated splitting of 2.2 meV is 
resolvable with the existing neutron scattering spectrom- 
eters. For completeness we note that the combined effect 
of the external magnetic field and the inter-layer coupling 

on the Goldstone modes is given by w - 



y 



- uJ 



V. CONCLUSIONS 

The excitations in the antiferromagnetic cores of su- 
perconducting vortices in the 50(5) model have been 
studied. By examining the existing literature on experi- 
mental results connecting to the values of the parameters 
of the model, we have found that the stiffnesses in the 
charge and spin sectors are nearly identical, p c ~ p s , and 
likewise for the susceptibilities, \c ~ Xs- This remarkable 
fact serves as good support of the idea of the existence 
of a 50(5) symmetry in the high T c cuprates. 

We have predicted within the 50(5) model that bound 
localized excitations exist when asymmetries arise be- 
tween the 7r sector of the parameters and the spin and 
charge sectors. If they exist, these excitations could be 
observed in neutron scattering experiments as side peaks 
to the already observed tt excitation, side peaks with an 
amplitude proportional to the number of vortices and 
thereby proportional to the applied external magnetic 
field. 

Finally, we have predicted the splitting of the degen- 
erate zero energy mode as a function of applied mag- 
netic field and the inter-layer coupling. The effect of the 
magnetic field is minute, only a few /J,eV, and thus not 
possible to detect with present day neutron scattering 
technology. The effect of the inter-layer coupling, on the 
other hand, is of the order of 2 meV and hence detectable 
in inelastic neutron scattering experiments. The test of 
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the existence of these core excitations would constitute 
a crucial test of the 50(5) model. The expected signal 
should only be present in the superconducting phase, and 
it should be proportional with the number of vortices, i.e. 
with the applied magnetic field. 
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